Indonesian Journal of Electrical Engineering and Computer Science 
Vol. 29, No. 2, February 2023, pp. 703~714 
ISSN: 2502-4752, DOI: 10.1159 1/ijeecs.v29.i2.pp703-714 0 703 


Overcome uncertainties of vertical take-off and landing aircraft 
based on optimal sliding mode control 


Wajdi Sadik Aboud', Hayder S. Abd Al-Amir’, Aseel A. Alhamdany’, Fahad M. Kadhim! 
‘Department of Prosthetics and Orthotics Engineering, College of Engineering, Al-Nahrain University, Baghdad, Iraq 
Department of Power Mechanics, Institute of Technology, Middle Technical University, Baghdad, Iraq 
3Department of Electromechanical Engineering, University of Technology, Baghdad, Iraq 


Article Info 


ABSTRACT 


Article history: 


Received Aug 7, 2022 
Revised Oct 11, 2022 
Accepted Oct 24, 2022 


Keywords: 


Particle swarm optimization 
Sliding mode control 
Stability 

Tracking 

Vertical take-off and landing 
aircraft 

Vibration 


This study aims to design a robust and optimal controller to overcome the 
problems related to the existence of disturbances and uncertainties during 
takeoff and landing operations of vertical take-off and landing (VTOL) 
aircraft. The dynamics are decomposed into two phase’s parts which are the 
minimum phase and the non-minimum phase. These two-part are controlled 
by proposing a robust nonlinear controller represented by sliding mode 
control (SMC). Also, the chattering effect due to the fast-switching surface 
in SMC is eliminated by utilizing a proposed sigmoid function which acts as 
the sigmoid function. The controller's main parameters are tuned optimally 
based on the particle swarm optimization (PSO) algorithm. In addition, the 
controller guaranteed the system stability based on the Lyapunov and Routh 
theories. The main output parameter responses represented by the 
positioning of the centre of mass and angle of rolling are determined with 
bounded control inputs. The performance of the proposed controller is tested 
by tracking VTOL parameters to the desired trajectories. The simulated 
results not only showed a significant tracking trajectory but also system 


stability guaranteed. In addition, the results showed an improved rate of 72% 
and 84% compared with those results obtained from the literature. 
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1. INTRODUCTION 

Fixed-wing combat aircraft is considered one of the most important tools of modern armies because 
of their combat, defence, and reconnaissance missions. The main important requirement for this type of 
aircraft is large and flat areas for take-off and landing. These areas are known as airports, whether they are on 
land or on ships. These airports are far from the battlefield and cannot always be provided near them. 
Therefore, there was a need to think about adding the vertical take-off and landing (VTOL) feature to fixed- 
wing aircraft as in helicopters. That supports these types of aircraft to land almost anywhere and eliminates 
the need for long airport runways. This helps to reduce the risks associated with the movement of the aircraft 
when taking off and landing, especially on battlefronts. The VTOL aircraft have many applications and 
usages such as reconnaissance, search and rescue, logistics and other fields. The most important feature is the 
use of secret places that can be hidden for take-off and landing. The landing gear of this type of aircraft is 
less complicated than that of fixed-wing aircraft. However, it needs flat land without ripples to land and takes 
off, which reduces its work in all types of terrain and limits its capabilities. It is noteworthy that there are 
numerous types of VTOL vehicles, which it was listed through a reference [1]. The most worrying of the 
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engineering designers is the necessity of having high stability of the plane for vertical take-off and landing. 
This is represented by preventing rotations around its own three axes (yaw-roll-pitch), especially in difficult 
weather conditions such as high wind. During the take-off and landing of VTOL aircraft, due to the presence 
of a vertical thrust downward, uneven forces appear on both ends of the wing in the vertical direction, which 
causes a roll moment reaction about the aircraft's longitudinal axis. The slight lateral acceleration of the 
aircraft will be generated as presented in the current interesting study. This phenomenon makes the plane 
unstable. The problem of the stability of the aircraft at vertical take-off and landing has been solved by 
introducing the error tracking system that occurs with the stability of the aircraft in the horizontal plane. 

To study the stability of the aircraft during take-off and landing, simplified two-dimensional 
mathematical models have been developed that simulate well the dynamic behaviour of the aircraft. These 
mathematical models have proven their efficiency and quality for various scientific applications over the past 
years. Then provide the available automated control algorithms for the desired results. Real experiments on 
control equipment are very expensive in terms of manufacturing technology as well as in terms of testing 
stages or trials. The presence of these models allows researchers to develop models of control units of various 
types to simulate the control of aircraft movement during take-off and landing and increase its stability [2]- 
[4]. These controller structures were built on one tracking control algorithm such as approximation input- 
output linearization [5], system inversion [6], linear optimal control [7], internal-model based approach [8], 
backstepping [9], and Lyapunov function technique [10]. To achieve the reality of dynamic systems, 
(certainty) coefficients are introduced to examine and develop the proposed controls to be more realistic 
before they are applied to the systems in real life. Researchers have introduced many control strategies during 
the take-off and landing of aircraft (VTOL), to make them more stable and safe under various conditions. 
The controller type proportional, integral and derivative (PID) has been used a lot with linear models because 
it is still really simple and easy to control and implement later as in [11]-[15]. Waslander et al. [16], a 
comparison of two nonlinear controllers based on integral sliding mode and reinforcement learning is 
presented. More recently, the studies focused on robust controllers when disturbances and uncertainties were 
available [17]-[21]. All control systems applied in the field of aircraft seek to make the angle of rolling close 
to zero and reduce the hover fluctuation of the aircraft during vertical take-off and landing as well as 
suppression of the undesired vibrations. Generally, there is a lot of work that has to be done using the 
available control theories to practical implementation and down to complex details. 

The essential motivation of this study is possessed from [16], [17], [20] because the stabilization 
issue of the VTOL at the centre of mass for the aircraft is still there waiting to be treated. This can be done by 
producing force control actions for the VTOL able to track the desired trajectories with a minimum error. 
And also able to eliminate unwanted bounded disturbances. Moreover, the characterization of the 
contribution of this research according to the comparisons with the literature that includes different types of 
nonlinear controllers for VTOL can be enumerated as follows: i) the control law for both the minimum phase 
and the non-minimum phase dynamic parts is developed based on very accurate analytical derivation using 
sliding mode control (SMC). In addition, the stability based on Lyapunov and Routh criteria is guaranteed 
and ensures optimal control forces for fast control action by using the particle swarm optimization algorithm 
technique; ii) very fast tracking trajectory and minimum tracking error for different initial states; ii) the 
chattering phenomena of the sliding mode control are eliminated by minimizing the output of the signum 
function with the minimum value as possible; and iv) the proposed controller is capable to overcome the 
effects of undesirable disturbances and parameters uncertainties effectively. 

In this study, the dynamic mathematical model of the VTOL based on Newton's law is derived and 
analysed. The derived dynamic decomposed into the minimum phase dynamic part and the non-minimum 
phase dynamic part. The first part is used to control the vertical dynamic of the flight and the latter part is 
used to control both the horizontal and roll dynamic of the flight. Optimal SMC with bounded input control is 
proposed to obtain trajectory tracking and stability. This work aims to provide a tracking automatic control 
with a high potential to control the movement of the aircraft during vertical take-off and landing. This was 
achieved by making the angle of rolling close to zero and reducing the hover fluctuation of the aircraft during 
vertical take-off and landing as well as suppressing the undesired vibrations. 


2. MATHEMATICAL MODEL OF VTOL 

According to Figure | and based on the principles of Newton’s law, and let the small coefficient, €9, 
which describes the relationship between the lateral force, €)!, and the rolling moment of the VTOL aircraft, 
the dynamic model is derived as [5], [6]: 
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Figure 1. Planar VTOL aircraft 


—mX = —-(T + 6, (t))sin(@) + & (1 + 62 (t))cos(6) 
—mY = —(T + 6; (t))cos(6) + €9 (I + 6 (t)) sin(@) — mg (1) 
1,6 = 1 + 6,(t) 


The dynamic model given in (1) can be rescaled with the following are define x = —X/g, y = —Y/g, u, = 
T/(mg), uz = l/Iz, € = Eglz/(mg), &(t) = 6;(t)/(mg), &2(t) = 62(t)/Iz. Then, the dynamic model can 
be written as (2). 


i= —(u, + &,(t))sin(@) +€ (uz + & (t))cos(@) 
je —(u, + &, (t))cos(0) + € (uz + & (t))sin(@) -g (2) 
6 =uz + &, (t) 


Suppose that x; = X,X%2 = X,y, = Y,Y2 = y. Then, the scaled mathematical dynamic equation of 
VTOL becomes as below; 


x=x 

X_ = Ge + & (t)) sin(@) + e(uz + & (t))cos(6) 

Y=Y2 

V2 = (u, + & (t))cos(@) + € (uz + & (t))sin(@) —g (3) 
0G=W 


@ =u, + & (ft) 


where x,(t) and y,(t) are represent the positioning of the centre of mass for the aircraft, 0(t) is the angle of 
rolling, X2(t) and y(t) are represent the linear and angular velocities respectively, the € is defined as the 
constant of coupling between the lateral force and the rolling moment and g is the gravitation acceleration. 
The control inputs, uw, and wz, represent the thrust and rotational moment respectively. The uncertainties of 
the systems are defined by 6), 62, and 63. It is clear from the model given in (1) that the aircraft model is not 
only underactuated but also in a non-minimum phase when €40 concerning the nominal output y. Consider, 


d,(t) = 6,(t) = —&,(t)sin@ + € &cosd 
dz(t) = 62(t) = §,(t)cos@ + e sind (4) 
d3(t) = 63(t) = & 


in this study, the differences between outputs of the VTOL system and the assumed desired trajectory, 
(%g,¥q), can be defined as the tracking error: 


é; =e 

é, = te + & (t))sin(@) + € (uz + & (t))cos(@) — Xq 

é3 = e4 (5) 
C4 = —(u + & (t))cos(@) +E (uz + & (t))sin(@) —g-Va 

d6=a 


@ =Uz + & (t) 
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from (5), the following expression can be driven [5]. 


iz i Bee ecos0 [i] a Bae aa de Xa (6) 
e4 cos@_ esin@1 U2 €,cos@ +é€&,cosé —2-Va 
The utilizing of the input state linearization is necessary because the relative degree for the vector of 


(5) corresponding to the output is [2 2]. This is can be done by adopting the control law presented in [22], 
[23]. So, 


[e2] = Bees eel atx (7) 
Uz] Lcos@ esinO! |A,+e+ Yq 


by the following steps, the 2, and A, can be designed [5]. 


—sin@ «cos0|~* —sinO cos 

| i =]1 Les (8) 
cos@_esind =c0s@ ~sin@ 

Therefore, converts (5) to, 

ey = e2 

é2 =A, + 6,(t) 

€3 = C4 

4 = Az + 52(t) (9) 

d=a 


@ =U, + 63 (t) 


then the following formula can be drawn from (7) and (8); uz = = A,cos0 +~ A,sind + = sin6 + 


~%qcos@ + - ¥qsin@. Based on that, the following equation is provided; 


é, =e 
és = 1, + 6,(t) 

e3 = C4 

d= 


O= - A,cosé + - A,sind + = sind + ~X4cos0 + ~¥jqsind + &(t) 
a new variable, 7, instead of w is considered for eliminating A, and A, in @ as (11) [22], [23]. 


X = Ew — e2cosO — e,sind (11) 


From (11), it can be seen that if e, > 0,e, > 0, then v > ew, and if v > 0, then w — 0, so that v 
can be used instead of w [22], [23]. Also, it can be seen that w includes A, and Ag, é, includes 2,, and é4 
includes A,, so by combining w, y, e, the v could be designed. Thus, the differential equation 7 can get with 
@, €2, and é,. Therefore, the elimination of 2,, Az is possible [22]. So, from (11), one can get (12). 


e=0= =(x + e,cos@ — e,sin@) (12) 


Then, 
XY = ew — €2c0S0 + e,OsinO — eysin@ — e,Ocosd 
=€ (Gar cosO + - A, sinO + 2 sind + - Xqcos 8 + - Yasin 0 + é2(t)) 
—(A, — &,(t)sin@d + e&,(t)cos@)cosé — (Az + &,(t)cosé) (13) 
+ €€,(t)sin9)sin@ + - (e2 sin? — e, cosO)(y+e2 cos@ + e, sin) 
= ~( €z sin@ — e, cos@)(y+e, cos@ + e, sin9) + gsin@ + Xqgcos 6 + yqgsin @ 


hence, the dynamics part of the non-minimum phase is obtained from (12) as (14). 
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v4 = q(e1, €2,€3,€4,9,X, Ya) (14) 


where ¥=[@ Y]", Ya = [Xa al, and, 


1 

Pie +e, cos@ + e, sin@) 

q(e1, €2, 3, €4,9,X, Yq) -_ 1 
at e,sin@ — e,cos@)(y + e,cos@ + e, sind) + gsin@ + Xg cos 6 + ygsin 8 


It is evident from the latter equations that the q(.) is related to [@1 2] at zero states but it does not 
to [€3 4]. Therefore, the following can be derived [22], [23]. 


Aq (€1,€2,€3,C4,9,XYa)|  _ 
ee Cr ee Tl 2x2 (15) 
4q (€1,€2,€3,€4,4,X% a) 
= ae ee ae Ilpx2 (16) 

Based on (15) and (16), it is clear that the model given in (10) can be decomposed into two parts. 
The first one is the minimum phase part that is associated with control of vertical flight dynamics and the 
second one is the non-minimum phase part that is associated with control of the coupled horizontal and roll 
flight dynamics. The two parts models are presented in (17) and (18) respectively. 


é3 = e,(t) 

ey = 5 

eo => 4 + 6,(t) 

é= “(x + e,cos@ — e,sin@) (18) 


X= - (e sin@ — e, cos0)(y + e2 cosO + e, sind) +(g+¥q) sind + Xgcos @ 


3. THE DESIGN OF THE OPTIMAL CONTROLLER 
This study uses feedback inversion represented by a conventional sliding mode controller to resolve 
the output-tracking problem for the minimum phase dynamics given by (17). The SMC is selected as (19): 


0, =cez+6é3, e3(t) =v - Ya,» €4(t) = V2 - Va (19) 


where c is a positive scalar parameter used to guarantee the system trajectory to hit the sliding surface when 
t — co, i.e. c > 0, then the controller is selected: 


Az = —ce, — h,sgn(o,) (20) 


where h, > |d, (t)|, and sgn is the signum function that is proposed in this study as a sigmoid function, 
sgn(o;) = c tanh(o;) where i = 1,2, to reduce or overcome the chattering effect related to the SMC [24]. 


The Lyapunov function is defined as V,; = ~o3, then the following can be derived [25], [26]. 


V, = 10, = (céz + é4) = 0, (ces te (A, +dz)) = 0, (ceg + (—ceg — hysgn(o,) + dz)) 
= 0;(—h,sgn(o,) + d2)) = c(—Ay lo; | + d0,) <0 (21) 


Therefore, for any differential output command yg , y(t) > vq , y2(t) > Yq as t > co. Thus the vertical 
dynamics and both the horizontal and roll dynamics of the aircraft do not affect each other because the 
minimum phase dynamics and the non-minimum phase dynamics are completely decoupled. Here, the main 
goal is to design SMC using the above sigmoid function to stabilize the system regarding the non-minimum 
phase dynamics. So, for (18), let denote, 


Hy = @2,U2=[e1 9 X) ervH my — My ep Se Ky (22) 
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then, (18) is written as; 


fy =A, +d, (0) 
: 23 
Hz = P(e1, €2, 3, €4, 9, X, Ya) ) 


e2 


where p(é1, €2, €3,€4,9,%, Ya) = Ly-+e2 cos 6-+e4 sin 6) 
Lurteg cos 6+e4 sin 8)(e2 sin 0—e4 cosO+Xq cos O+(Vqtg )sind 


based on Taylor's expansion, (23) can be written as, 


_ — 9P op T 
Bo gel 2 ore, oat [er @ Xa)" + M(e1,€2, 3, €4, 4, X, Ya) (24) 


= Assy Uy + Ass U2 + H(e,, €2,€3,€4, 8,%, Ya) 


where, 
0 0 O 
g 
M(e, €2,€3,€4, 8%, Ya) = p(er, €2,€3, 4, 0, xX Ya) mn Ass U2 - Avsa My 
The sliding variable for (23) is defined as, 
(25) 


02 = ty, — Mun 


where M=[M, ™ 2 Ms] and M is selected in such a way that let A >, + A 2,M Hurwitz. According to 
(25), the sliding time, ts, exists if and only if the sliding mode exists. Therefore, for t => ts which let 2 = 


O, then the following can be derived based on (24). 


fly = A ny Miz + A 99 My + M(€y, €p, €3, C4, 9X, Ya) (26) 
- (A 2M+A 22) Hz + II(e€1, €2, €3, €4,9,X, Ya) 


The controller parameters are selected as: 


A, = Mp (ey, €2,€3, €4,9,X, Ya) — hz sgn(o2) (27) 
where h, > |d,(t)|. Again, the conventional Lyapunov function is V = <o8, Therefore, 
V = 02 (tt, — Min) = 02(A, + d,(t) — M tz) 
= 021Mp (ey, €2, €3, €4, 9, X, Ya) — he sgn(oz) + d,(t) — Mp(ey, €2, €3, €4,9,X, Ya)} (28) 


= dn {d,(t) — hy sgn(o2) } = 02d,(t) — hy |o.| < 0 
now by executing the vector M which is equal to M = [m, mz m3] and let A> + A 21M to be Hurwitz 


Then, the following get. 


0 0 O 1 0 0 O mM, M2 M3 
Ay2+A2iM = f O e*)+fet{[m, mz m3] = [ O e*l+letm, etm, etm, 
0 g O 0 0 g O 0 0 0 
mM, m2 m3 
= c tm, elm, e} rem (29) 
0 g 0 


Hence, 
|sI — (42. +A 21M)| = 53 — (m, + €-1m,)s? — g(e71 + €-1m,)s + gem, (30) 


based on the rule of Routh, the following relations are obtained (31). 
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elm 


Mm, > 0,mz < —Em,, m3 < — (31) 


m,t+e-1mz 


The SMC control parameters (c,M,and h; where i = 1,2 ) are selected optimally by using the 
particle swarm optimization (PSO) algorithm technique [27]. The idea behind the PSO technique is that the 
best solution can be obtained by simulating the movement of the birds. This technique adopted the particles 
as a population of the individual which flies over the space of solution finding the optimal result or solution. 
The particles are moving around the search space and each one of these particles has its position and velocity. 
Hence the evaluation of particles is based on the closest one to the optimal solution and that can be done 
according to the function named fitness function [27]. 

The PSO algorithm depended on two main values. the first one is the previous best value named 
pbest which is related to the particular particle. The second value named gbest is the best compared to all 
the particles pbest in the overall swarm. The SMC with four weight parameters is arranged in an array to 
represent the particles. Then initialized all particles randomly and the updating followed afterwards based on 
(32)-(35) [28] to obtain the optimal values for the sliding mode controller parameters in all cases that 
encounter the work of the aircraft. 


Ac;** = Ac} + kyrand,(pbest; — c;) + kyrand,(gbest' — c}) 


been aet ; (32) 
cj*t = Ac; + Acj** 
AM;** = AM} + k,rand,(pbest} — M;) + k,rand,(gbest' — Mi) 3 
Ah ae = Ah + k, rand, (pbest} —h i) + k,rand2(gbest! —h A) 
hist = Ahi, + Ah it a 
Ah 1 = Ah}, + kyrand,(pbest! —h},,) + k,rand,(gbest' —h$,) 
2j 2j 1 1\P J 2j 2 29 2j (35) 


AS = AhG, + Ah Gh" 
where j = 1,2,3...Np, Np defines as the number of particles, cj, Mj, hi i and hi, ; define as the particle weight 
j at i iteration, k; and kz define the constants acceleration which has a value of 2. The values rand1 and rand2 
represent the random numbers where used here between 0 and 1, pbest' defines as the best previous weight 
of j" particle, and gbest is the best particle among all the particles in the population. 


The number of dimensions in PSO is equal to four due to the four main controller parameters of the 
SMC. The model estimation is based on the criterion of the mean square error function as (36), (37). 


1 GN. 
E, = nym (ext + 


36 

Ey = = YnP(eo(i + 1)")? we 
Np 

E3 = Sy? (e3(i + 1)")? 

Pane , (37) 

E,y= xpum (eat Ly") 


The values of the mean square error (E, — E4) must be less than 1e~° to satisfy the optimal tuning. 
Table 1 presents the parameter values adopted for the solution of the PSO of the current study. 


Table 1. PSO parameters 


Parameter Value 
k; 1.3 
kp 1.3 
No. of particles 4 
Population size 100 
No. of iteration 66 
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4. SIMULATION RESULTS AND DISCUSSIONS 
In this section, the analytical derivation of the system given in (3) and the proposed SMC that was 


derived in (7), (20), and (27) are verified by using the MATLAB platform. The first step is that the main 
parameters of SMC are selected by try and error procedure. The second step is that the main SMC parameters 
are tuned optimally based on the swarm optimization technique available in MATLAB options. The 


following optimal parameters are selected. 
= = 7,.=[° = 
M=[5 -14 —6], c=50, hj = Gal where i = 1,2 


The coefficient of coupling ¢ is taken in the current study as 0.5 where this value produces a non- 
minimum phase system of VTOL aircraft that is strong [29]. The trajectory of the desired output is assumed 
as Xq = Rsin(wt), yg = Rsin(wt), where R is the amplitude and o is the frequency of the desired input. 
These values, R and w, are taken from 0.1 to 0.2 and | to 1.2 respectively to simulate a case of obstacle 
avoidance [29]. The initial conditions are selected as x(0) = [0.5 0.5 0.5 — 0.05 1 0.5]". 

To represent the present robustness of the proposed controller, the state tracker and generator are 
assumed as €, = 0.5sint and €, = 0.5cost and performed the simulation with these values. Notice that &, 
and €, are inertial parameters, which are not important to know their values practically [6]. Figure 2 
illustrates the block diagrams of the proposed controller scheme in MATLAB Simulink. 

Figures 3-5 show the output tracking simulation results. Figure 3 presents the tracking trajectory of 
the positioning and the velocity of direction x to the proposed desired trajectory. It is evident that the time 
required to have zero error is around 2.8 s. Figure 4 presents the tracking trajectory of the positioning and the 
velocity of direction y to the proposed desired trajectory. The time required to have zero error is the same as 
that for x-direction. Figure 5 shows the tracking trajectory of the angular positioning and the angular velocity 
to zero, which is the proposed reference trajectory. It is obvious that the time required to have zero error of 
the angular positioning and the angular velocity is around 2.4s and 2s respectively. 
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Figure 2. Simulink representation of proposed optimal SMC of VTOL system 
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Figure 3. The positioning and velocity responses in x-direction versus the desired responses 
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Figure 4. The positioning and velocity responses in the y-direction versus the desired responses 


6/ 

i=) 

i 

=| 

S \ imulate 

=) — reference 0) || 

4 ae — \ 1 1 
0 5 10 15 20 25 30 35 40 45 50 
time(s) 

3 

3 | = 

= —————EE 

§ -1) [ly se simulate w 

S i — ref w 

© \ } reference 

5 204 /— i i i i i i J 
o~ 5 10 15 20 25 30 35 40 45 50 


time(s) 


Figure 5. The roll angle versus the zero responses 


To monitor the stability of the closed-loop feedback nonlinear optimal SMC for the VTOL system, 
Figure 6 presents the phase-plane plot for concentrating upon the chattering influence which is discarded due 
to these causes. The signum function utilized the Sigmoid function, sgn(o;) = c tanh(o;), in this study and 
the control parameters values are calculated by adopting the Particle Swarm optimization algorithm. Figure 7 
shows the behaviour of the sliding variable (surface). It is evident that the effect of chattering is highly 
eliminated at the fast-switching surface. This is due to the use of the signum function which leads to reducing 
the amplitude of the output function then the responses of VTOL approached the desired responses at the 
inverted position. It can be noticed that the control inputs, uw, and uz, presented in Figure 8 are rather 
aggressive behaviours that take place in relatively short periods as well as they are bound. In addition, the 
selected controller and the controller behaviour are presented in Figures 9 and 10 respectively. 
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Figure 9. The selected controller behaviour Figure 10. The controller behaviour 


For validation purposes, the essential results of this study were compared to some of the literature. 
Despite the different techniques and the aim in [29], [30], the responses are compared for both positions and 
velocities. The results for configuration tracking and velocities are shown in Figure 11 according to [30] 
while the result of roll angle variation is illustrated in Figure 12 according to [29]. By comparing these results 
with results presented in Figures 3-5, it is evident that the time required for steady-state is improved by 74% 
and 71% for configuration and velocities respectively according to [30]. In addition, the comparison of roll 
angle response is 81% less in Figure 3 compared to the results of [29] as shown in Figure 12. These 
improvement rates show the effectiveness of the proposed controller in this study. 
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Figure 11. The actual configuration responses versus the Figure 12. The actual roll angle response 
desired responses [30] versus the desired response [29] 
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5. CONCLUSION 

In this study, the dynamics of VTO are analyzed based on Newton’s law. The SMC is proposed to 
track the desired trajectory. The SMC uses a proposed sigmoid function to eliminate the chattering effect. 
Also, the controller is tuned optimally by using the PSO algorithm technique. In the designed optimal SMC, 
not only the tracking trajectory is obtained but also the stability is guaranteed based on Lypinove and Routh's 
theories. The MATLAB Simulink is used to execute the proposed controller. The results show the significant 
output responses for all main parameters assumed. The error goes to zero in a very significant period not 
reaching 3 sec. The validation results showed an improvement average rate reaching about 72% and 84% for 
reduction into zero steady-state responses according respectively. Not only have these significant responses 
improved but also stability is achieved when the disturbances are available. Also, the chattering which is an 
SMC disadvantage property is highly eliminated as well. 
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